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We revisit the hard-spheres lattice gas model in the spherical approximation proposed by 
Lebowitz and Percus (J. L. Lebowitz, J. K. Percus, Phys. Rev. 144 (1966) 251). Although no 
disorder is present in the model, we find that the short-range dynamical restrictions in the model 
induce glassy behavior. We examine the off-equilibrium Langevin dynamics of this model and study 
the relaxation of the density as well as the correlation, response and overlap two-time functions. 
We find that the relaxation proceeds in two steps as well as absence of anomaly in the response 
function. By studying the violation of the fluctuation-dissipation ratio we conclude that the glassy 
scenario of this model corresponds to the dynamics of domain growth in phase ordering kinetics. 

I. INTRODUCTION 

The nature of slow relaxation in frustrated systems has seen a large increase of activity-in the recent years. In 
particular much effort has been done in the study of relaxational dynamics in spin glasses^. These are disordered 
systems where slow relaxation appears as a consequence of the frustration induced by the disorder. At low enough 
temperatures the system explores a rugged free energy landscape with high energy barriers, hence dynamics is slowed 
down. But nature offers a large variety of systems where dynamics can be exceedingly slow in the absence of disorder. 
Structural glasses are examples in this class of systems. When fast cooled, glasses do not crystallize and leave the 
undercooled liquid line when the typical relaxation time exceeds the inverse of the cooling rate. A possible scenario for 
the origin of the glassy statain the absence of quenched disorder has been recently proposed in the framework of some 
solvable mean-field modelsou. These models show the existence of three characteristic temperatures: the melting 
first-order transition temperature T m , a dynamical transition temperature Td reminiscent of a spinodal instability 
and a thermodynamic glass transition T s where the configurational entropy vanishes and replica symmetry breaksO. 
A detailed study of the dynamical equations of these mean-field modelsD shows that these indeed correspond to the 
mode coupling equations of Gotze and collaborators for the structural glass problemEl. One of the main results in this 
approximation is the existence of a dynamical singularity (T<j) where dynamics is arrested and ergodicity broken. The 
existence of this singularity relies on the mean-field character of such approximation and it would be highly desirable 
to understand how to include short range effects in the theory in a systematic way. 

To address this problem it is useful the study of solvable microscopic models for a better understanding of the 
mechanisms responsible for slow relaxation, the universal properties of the off-equilibrium dynamicsE3£3 as well as 
its connections with the Mode Coupling approach. In this paper.. we study the off-equilibrium behavior of a solvable 
spherical model introduced long time ago by Lebowitz and PercuaJ. The model has not disorder built in and is suitable 
to study the role of short-range dynamical constraints in glassy systems. We will show that the off-equilibrium behavior 
of this model shares the same common features (and in some sense, belongs to the same universality class) as those 
mean-field spin glass models characterized by the absence of anomaly in the response function. Quite remarkably the 
glassy scenario of the model corresponds to the domain growtbEll. We should stress also that the model introduced 
by Lebowitz and Percus is not enough to provide a complete description of the glass transition where phenomena like 
stretching and the characteristic two step form above T c in the liquid (disordered) phase are absent. IjUjhc language 
of Mode Coupling theory this is due to the absence of discontinuity at T c of the ergodicity parameter&Ej. 

But still this model deserves its interest. The new feature which this model furnishes to the study of glassy dynamics 
is the effect of dimensionality and short range dynamical constraints in the mechanisms foe. glassy dynamics. This 
kinetic effect has been proposed in the past as a possible explanation of the glass transition^. The Lebowitz-Percus 
model (hereafter referred as LP model) incorporates the effect of spatial correlations in the system. This allows for 
an extension of the dynamical equations to include the wave vector dependency. 

The paper is divided as follows. In the section II we define the model as well as its thermodynamic properties and 
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write down the dynamical equations in the Langevin approach. In section III we close the dynamical equations for 
one time quantities. Section IV presents the solution of the dynamical equations for the correlation, response and 
overlap two-time functions. Section V discusses some numerical experiments of the model which allow for a clear 
identification of the glassy scenario. Section VI discusses the nature of relaxation processes as well as violations of 
the fluctuation-dissipation relation. Finally we present the main conclusions. 



II. THE LEBOWITZ-PERCUS (LP) MODEL 



Let us consider a lattice gas in the spherical approach. Lattice gas models arc defined on a lattice of finite dimensions 
where local densities p(x) are attached to every site of the lattice. In the hard spheres lattice gas every site of the 
lattice can he either occupied or empty. In this case the local densities p(x) can take only the values or 1. Lebowitz 
and Percuso relax these conditions on the local densities allowing them to take any continuous value but satisfying a 
global spherical constraint. Following their work we also introduce the additional restriction that the density-density 
correlation function between nearest neighbors vanishes. This restriction mimics some kind of extended hard core. 
So the restrictions on the system are, 

a 1 ^^p(x) 2 -^p(x)=0 (1) 

X X 

<7 2 = ^^p(x)p(x + q) = (2) 

x q 

where x is a discrete variable that runs over the sites of a D-dimensional lattice and q are the vectors that join a 
site to its nearest neighbors. The thermodynamics of this system is very well known for the general case considered 
by Lebowitz and PercusQ where a pairwise interaction potential was also studied. For simplicity, here we will not 
introduce any interaction potential between particles or additional restrictions over the local densities. We will see 
that in this case the dynamics of the model is simply enough to be solvable displaying a non trivial relaxational 
dynamics. Short range dynamical constraints are such that relaxation turns out to be slow when a reorganization of 
local densities is necessary to reach the equilibrium density of the system. We start by solving the thermodynamics. 
After, we define the dynamical equations of the model and find the stationary solutions of the dynamics. 



A. The thermodynamics. 

In the grand canonical ensemble the partition function of the LP model is given by, 

oo 

£ GC = $>"Q„(/3) (3) 

n=l 

and Qn(P) is the partition function computed in the canonical ensemble, 

Q n (f3) = £ ^ 1)0 <W (4) 
{«}' 

where u\ and o-i have been defined previously in eqs.([l],^|). The prime in the sum of eq.(||) indicates that the set 
of densities {p{\ satisfy the global constraint Y2i=i Pi — n where n is the total number of particles and N the total 
number of sites in the lattice. Note that the thermodynamics of this system is determined solely by the entropy 
because there is no interaction energy between the particles (like in the general case of hard spheres systemsE3). 
Substituting the canonical partition function eq.(|^) into the grand partition function eq.(^) we obtain in the large 
n, N limit with n/N fixed, 

oo N 
Pl=— oo 2 = 1 

Introducing the integral representation of the delta function we obtain, 
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2 = n df>l e.xp[f3^C{x u x v )pi{l - p v ) + ^2ipi] (6) 
J i=\ i.v i 

where z is the fugacity z — exp(7) = exp(/x/3) and C{xi, Xi>) = X(xi — xi>). In our case, X(xi — xy) = Ao when xi = xy 
and X(xi — Xi>) = Ai when xi — xy — q, otherwise X(xi — xp) = 0. 

The quadratic exponent can be diagonalized using the Fourier transform and the partition function can be exactly 
evaluated in terms of the Lagrange multipliers. The equations for the multipliers can be obtained noting that 
= cf\ and dl g^ Z — (T 2 which yield eq.(|l|,||). Using these constraints we obtain the following equations for the 
Lagrange multipliers, 

1 _ T cosjkx) = 

AT y 2A ( k ) 

A(k) = A + 2Ai^cos(k.q) (8) 

1 

where x = or x = with i = 1, D. In what follows we will drop out the brackets in the notation for the average 
density writing p instead of < p >. These equations can be readily solved. In the disordered high temperature phase 
eq.(Q) takes the simple form, 

1 f 2n D Tcos(k_x) 

Wf 5 Jo 2X = P{S(X) P) (9) 

(10) 

where again x = or x = q^, (1 < i < D). A condensation phase transition is found above two dimensions. This result 
can be easily inferred after examination of eqs. (Q),®. The phase transition corresponds to the condensation in the 
K c = 7r(l, 1, .., 1) direction due to the positivity of Ao and Ai. At the transition point A(K) = 0, i.e. Ao = 2Ai-D = A. 
The critical temperature T c and the value of the Lagrange multiplier A are solutions to the equations, 

1 f 2 * n d T c cos(k-x) ... . , 

' rfk nw, i 1 ^ ft \T = P(^( X ) _ P) (11) 



(2n)D J 2A(l + -LE 9 cos(k.q)) 

x = or x = qi, (1 < i < D). Below the critical temperature, condensation on the mode K = 7r(l, 1, .., 1) develops and 
the Lagrange multipliers remain fixed to their critical values. The condensed mass ax as a function of temperature 
satisfies the equation, 



a K cos(K ■ x) + —- jr / dk D — '— — = p(5(x) -p) (12 



L_ r, h D rcos(k-x) 
(27r)^y 2A(l + iE,cos(k.q)) 

(13) 

where x = 0, q^, (1 < i < D). Thejaechanism for the condensation transition in the LP model is the same as in the 
spherical model of Berlin and Kadi] and appears only for D > Di = 2 which is the lower critical dimension of the 
model. For D — 3 previous equations can be solved and we find T c = 20.158. We should also note that the entropy 
of this model diverges as log(T) at low temperatures, hence violates the third law of thermodynamics. This implies a 
finite specific heat and a finite compressibility at zero temperature. This is not a serious drawback since this is related 
to thej-cpntinuous character of the local densities. To suppress this undesirable effect one should include quantum 
effcctsM but this is beyond the scope of the present paper. 

B. The dynamics. 

The effective Hamiltonian of the LP model which takes into account the restrictions (|l|) and (||) imposed on the 
system reads, 

ffe//(W) = Ao (*) 5>« 2 - P(x)) + A x (t) J2 E PMPfr + «») ( 14 ) 

X X l — l 
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where n v is the number of nearest neighbors, and qi are the vectors that join each point with their nearest neighbors. 

We propose the following differential equation for the dynamical evolution of the density in an open system that 
can interchange particles with a reservoir, 

f!^,-*f« t ,^l (15) 

where (3 is the inverse temperature and p is the chemical potential of the thermal bath. ?y(x, t) is a white noise 
uncorrelated in time and space such that (?y(x, t)) = and (f](x, t)r](x.' , t')) = 2 TS(x — x.')5(t — t') where the brackets 
indicate realizations over the thermal noise. From now on, we will only specify the temporal dependence when two 
different times appear dropping the explicit time dependence for one time dynamical quantities. 

We stress that in the LP model there is no energy but only entropy and that the role of the parameters Ao(i), A x (t) 
in the effective Hamiltonian eq.(|l4j) is to make the dynamical evolution for the densities to fulfill the dynamical 
constraints eqs(|l|,^|) for all times. Note that the role of the chemical potential (p > 0) in the dynamical equation ( |l5| ) 
is to increase the local densities in the lattice as much as possible. Obviously the density cannot increase indefinitely 
because of the dynamical constraints. Starting from an empty lattice the dynamics turns out to be slow when the 
density of the system approaches the equilibrium density. This slowing down is a purely entropic effect and is direct 
consequence of the decrease in the number of available configurations in pha se space. In the rest of the paper and 
without lost of generality, we can set the chemical potential equal to 1 in (|l5|). The dynamical equations read, 

p(x) = 1 - A (2p(x) - 1) - A x ^(p(x + qj) + p(x - qi )) + r?(x) (16) 

i 

Due to the spatial translational invariance, it is easy to diagonalize the system using the Fourier transform, 



p(k) = -Lj^exp (ik.r)p(r) 



(17) 

r 

The Fourier transformed global restrictions (jlj) and (§) are, 

( 7 1 =^|p(k)| 2 -ViVp(0) = (18) 

k 

a 2 =^^|p(k)| 2 7 (fc) = (19) 

k 

with 7(k) = jj J2h=i cos (k.qi). in terms of the Fourier components the effective Hamiltonian eq.(|l4|) is diagonal, 

H(p(k)) = Ao(t)(£ IWI 2 - VNp(0)) + 2D Ai(t) £ |p(k)| 2 7 (k) (20) 

k k 

and the equations of motion for the Fourier transformed densities are 

p(k) = VNS(k) - A (2p(k) - VnSQl)) - 4\iD p(k) 7 (k) + ry(k) (21) 

*K k ) = -L5> x P( ik - x M*) (22) 

* X 

This set of equations involve different uncoupled Fourier modes and can be formally integrated for each mode like 
in the Sherrington-Kirkpatrick spherical modelt3. Here we will follow a different strategy and write a linear partial 
differential equation for the evolution of the density. This is done in the next section. Before we will find the stationary 
states of the dynamics. Due to the absence of disorder we can also investigate the existence of crystal states in the 
system. 



1. Stationary states 

The stationary states of the model can be obtained by setting to zero the time derivative in (|2l]) and multiplying 
the equation by the complex conjugate of p. This yields, 
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= T 2 + A " + k) 4 nf°,:," <»> 

2Aq + AD Ai7(k) 



where we have used the regularization of the response function^, 

')p(t) >= 2T6 

< v(t'W) >= T . (25) 



lim < r]{t')p{t) >= 2TQ(t - t') (24) 



Using that |/5(0)| 2 = Np 2 , we note that 



2A + 4DAi y N 
The restrictions, o\ = and er 2 = eqs.(lqJl9|) yield 



P= +0(4)- (26) 



iV^2A +4A 1J D 7 (k) 
T 7 (k) 



iV^2A + 4 J DAi 7 (k) ^ 1 ' 

These are the equations originally derived by Lebowitz and Pcrcusi for the thermodynamics described in section 
II. A. At finite temperature we find that the only stationary solutions are given by the equilibrium states. This is not 
true at zero temperature where ergodicity is broken. In this case we expect the appearance of several crystalline states 
which nevertheless are metastable at finite temperature. Note that the phase transition in this model is different from 
the usual structural glass transformation where there is a melting first-order phase transition from a liquid to a crystal 
state. In the LP model the transition is a condensation one for D > 2 without any latent heat. 



2. Crystalline states 

At zero temperature, we note that there are many stationary states (its number being proportional to the size of 
the system N) . These can be obtained by setting T — in eq. ( |2l| ) which yields, after multiplying the equation by the 
complex conjugate of p, 

(2A + 4 J DA l7 (k))|p(k)| 2 =Np6(K)(l + \ ) (29) 

For k = this equation yields (2Ao + A\\D)p = (1 + Ao) but for k ^ we find different solutions depending on the 
value of k and — k where the term (2Ao + 4AiD 7 (k)) vanishes. Then the crystal states are characterized by a non 
vanishing value of p(0), p(k) and p(— k) = p(k)* , all the other modes being zero. The first term is the average density 
of particles while the second and third ones yield the crystal configuration. We still have to impose the conditions 
a% = and 02 — eqs. (|l^, [l9]). Then we obtain that 7 (k) should be smaller than zero for such a solution to exist. 
We also obtain that Ao = 1, and 2 Ai D — ~^jy (with k different from zero) and the equilibrium density is given by, 

7(k) rtm 
P = ~7T\ — 7 30) 
7 (k) - 1 

The number of stationary states is then proportional to the volume of the lattice since for each value of k such that 
7(k) < there is a stationary state. We will not extend further on details about the crystalline ground states but 
limit our discussion to the ID case. In this case, the simplest way to construct crystalline states is to assign a density 
equal to one to each point in the lattice every p sites. If p is a prime number, in the resulting periodic configuration 
only p(k — 0), p(k = — ) and p(k = 2tv(1 — -)) are different from zero. Then, the condition 7 (fc) < implies that k 

lies in the interval (§, ^). So, only the states with p = 2 (maximum filling of the lattice p = 1/2) and p = 3 (partial 
filling of the lattice with p = 1/3) are crystal states and fulfill the dynamical restrictions. All other periodic states 
with a larger value of p have j(k) > and decay to a new fc-state with j(k) < 0. 
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III. DYNAMICAL SOLUTION FOR ONE TIME QUANTITIES. 



To study the relaxation towards the equilibrium, we will focus our attention in the one times quantities as the 
density. We will need them later to study the evolution of two time quantities (correlation function, response function 
and the overlap among replicas). Dynamics is such that the global constraints a\ = (jl^) and a 2 = ( [l9| ) are 
satisfied for all times, hence di — and <72 = 0. These conditions determine the time dependent value of the Lagrange 
multipliers. 

2p-Ao + 4Z)Aip-l + 2i? = (31) 
X p - 4 Ai D T 2 + p + 2 i?i = (32) 

where p is the average density (p =< p >). The quantities T n are correlations density-density, and the R n quantities 
are correlations density-noise. They are defined as follows, 



N 

k 

= 4E(?W(k) +rXk)p*(k)) 7 (k)" (34) 

k 

where n is any integer number. Note that these quantities are invariant under translations which is the main symmetry 
of the effective Hamiltonian eq.(|l4|). We will see the usefulness of all these quantities later on. We can do some remarks 
on the values and physical meaning of some of them. First of all, we note that T is (due to the spherical restriction) 
equal to the average density p. T\ is proportional to the first neighbor correlation, which is equal to zero. T2 is some 
kind of second neighbor correlation. This is a time dependent quantity that needs to be calculated in order to solve 
the evolution of the density. We will see that the time evolution of the quantity T2 depends on T3, that T3 depends 
on T4 and so on. In the thermodynamic limit and using the regularization of the noise-field correlation (|2^) we find 
for the R n , 

k 

which are time independent quantities vanishing for odd values of n. 

We are now interested in the evolution of the average density of the system. Taking into account that (77) is of 
order O(^); we get, 

p= 1 - A (2p- 1) -4D\ lP (36) 

In order to solve the equations for the evolution of the density (|3l],|3^,^36]) , we need to know the dependence of T2 
on time. It is easy to write the dynamical equations for the T„, 

f n = 2p(l + A ) - 4A T„ - 8 D \ x T n+1 + 2R n (37) 

As previously said, each T n depends on T„+i. To close this hierarchy of equations, we multiply all of them by 

±x n and sum over n. Defining g T (x, t) = J2^=o Tr 7 "™ W and usin S tnat ^aT^ = X^=o TT^+iW' we S et a partial 
differential equation for gr(x, t), 



= _ g D Ai agrg& _ ^ gT ^ t) + 2p {l + \ Q ) e * + 2g N (x) 

9n (x) = J2 n ij xnR n = 7? Ek ex P ( x 7( k )) 
We can formally integrate this partial differential equation using the method of the characteristic curves, 



(38) 



gr(t, x) = g (x - /* 8 DXidt') exp (- f Q 4\ dt')+ , . 

/ *[2p(l + A ) exp (x - J* 8 D Aidt") + 2g N (x - f* 8 D Ai*")] exp (- J* 4X dt")dt' 1 ' 



where the go function is determined by the initial condition, gx(0,x) = go(x). Once the value of gx(t,x) is known, 

l g T {x,t) 1 



we can obtain the different elements of the hierarchy by taking derivatives with respect to x, T n (x, t) — 



G 



IV. THE HIERARCHY OF EQUATIONS FOR TWO TIME QUANTITIES. 



In this section, we write the dynamical equations for the correlation, response and overlap functions. These are 
defined by, 



G(t, t', x) = 1 <°( x ' + x > tMx'-.t') = ^ ]T /5(k, t)f}(-k, if) exp (zk.x) (40) 

x' k 

C(f, t', x) = l^ p(x' + x, i)p(x', = 1 £ p(k, t)p(-k, t') exp (ik.x) (41) 

x' k 

= ^E'° 1 ( x '+ X '*)' 92 ( x '' t ) = ]^X!piCM)P2(-M) (42) 

x' k 

C(t, t') (^) is the density-density correlation function and measures how fast the configurations decorrelate in time. 
The response function ( |40|) measures the change of the local density in a point of the, lattice at time t when the chemical 
potential is locally changed in another point of the lattice at distance x at time niJ. Finally the overlap function (|42|) 
measures how fast two different copies of the system (initially in the same configuration at t', i.e. pi(x, t) = P2{x, t')) 
decorrelate in time when submitted to different realizations of the noise for t > i1l3o. From now on and in the rest 
of the paper we will take the convention t > t' . To obtain their time evolution we proceed as in the case of the density 
and define hierarchies for two time quantities as follows, 

G„(t, t', x) = 1 £ p(k, t)fj(-k, t')j(k) n cos (k.x) (43) 

k 

C n (t, t', x) = 1 P(k, t)p(-k, t')7(k) n cos (k.x) (44) 

k 

Qn(t, if) = i Yl P^ k > *)^(-k, t)7(k) n (45) 

k 

Using ( pi] ) we get the following equations, 

BC 

-^(t, t', x) = -AD X 1 (t)G n+1 (t, t', x) - 2X (t)G n (t, t', x) (46) 

BC 

-^{t, t', x) = -AD Ax (t)C n+1 (t, t', x) - 2A (t)C„(* ) t', x) + p(i')(l + A (i)) (47) 
dO 

-j^(t, t') = -8D Xx(t)Q n+1 {t, t') - AX (t)Q n {t, t') + 2p(t)(l + A (i)) (48) 
For the response function we have < p{t')rj(i) >— 0. Defining the following generating functions, 



r G (t,t',x,y) = V^-G n (t,t , ,x) (49) 
*■ — ' nl 

n=0 
oo n 

r c {t,t',x,y) = V^-C„(t,t',x) (50) 

oo n 

rQ(t,t',y) = ^^ T Qn(t,0 (5i) 



nl 

n=0 



we can close the hierarchies (|49|),(|50|) and (|5 

-^(t, t', x, y) = -ADX x (t) (t, x, y) - 2A (t)r G (f , i', x, y) (52) 

-^L(t,t',x,y) = -ADX x {t)-^-{t,if,*,y) - 2X (t)T c (t, if, x, y) + p(t')(l + A (i))e y (53) 

-^(t, i', y) = —SD A x (*)-^(t, 2/) - 4Ao(t)r Q (t, f , y) + 2p(i)(l + Ao(t))e» (54) 
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The initial condition in equations (||||,[44],[l5|) is obtained setting t = t' . Using the regularization of the field-noise 
correlation at t = t' (p5|), the initial condition for the response function Tq^ (i', x, y) is given by 

(if, x, y) = r G (t', i',x,y) = E^^E ^ k ) n cos ( k - x ) ( 55 ) 

n=0 ' k 

For the correlation function, we have, 

00 „ 

(f , x, y) = T c (t', i',x,y)^^-^ \p(k, t')| 2 T(k) n cos (k.x) (56) 

n=0 ' k 

In particular, for x = we have (i',0,y) = gT{t',y). So we must use the generating function for the density as 
the initial condition for the generating function of the correlation. 

The same initial condition needs to be used for the generating function for the overlap between replicas, 

T$\t',y)=T Q (t',t',y)=g T (t\y) (57) 

Equations ( p2|j5^j5^ ) with their initial conditions ( p5| , p6| , ^7| ) can be formally solved as we did for the generating 
function for the hierarchy of the density ( |38"1 ) . 



3. The equilibrium solution. 

Using the integral expressions for Tc, Tq and Tq, we can find the equilibrium values for the two time quantities. 
In equilibrium, p, Ao and Ai are time independent, and the i nteg ral expressions can be simplified. Imposing also the 
equilibrium solutions as initial condition and using equation (|23|) we get the following results, 

reqf . s 1 ^ Tcos(k.x)exp(( 2/ -4JA 1 (t-^)) 7 (k)-2Aoft-t')) 2 ^ 

T c (t,t,x,y) = -^ 2A + 4£A l7 (k) +P 6Xp(2/) (58) 

T e J(t, t',x,y) = jjJ2 T cos ( k - x ) ex P ((» - 4 D Al (* - OfrOO - 2X o(t- 0) (59) 

k 

T e Ht ^^-iy T cos (k.x) exp ((y — 8D X, (t - Q) 7 (k) - 4 Ap(t - Q) , 

l Q (t,t,x,yj- 2A + 4£A l7 (k) + P exp (y) (bU) 

It can be readily seen that the equilibrium solution is time translational invariant and that the Fluctuation Dissi- 
pation Theorem (hereafter denoted as FDT) for the generating function is satisfied, i.e. arc( -g* / ' x '^ = Tc(t, t' , x, y). 
This implies the validity of the FDT for any of the elements of the hierarchy of equations. In particular, for the usual 
response and correlation functions we get dC °Q^ = Go(t,t',x). Also we note that T^(2(t—t'), x, y) — Tq (t— t' , x, y) 
as expectec@. 

In the thermodynamic limit the previous discrete sums in eqs.(|5^j5^,^) become integrals. In the condensed phase 
we obtain the following expressions, 

Vq (t — t', x, y) — ok cos(K • x) exp (y) + 

1 f 27r Wlf D T cos(k-x) exp ((y-4 D \ y (t-t')) 7 (k)-2 Ap (t-t 1 )) o pYri /, A 

ak 2A (i + i^ gCOs(k . q) ) P ex PW 

r£?(* - f, x, y) = a K cos(K • x) exp (y)+ 
^ J 2 " dk c T cos(k • x) exp ((y - 4 D X, (t - f')h(k) - 2 A (t - f )) ( } 

r Q (* ~t',y) = a K exp(y)+ 

1 f 27r A D T exp ((y-8 fl A t (t-t')) 7 (k)-4 A (t-t')) _ 2 PYT1 /,A 
P^p-JO flK 2A(l+i£ q CO S (k.q)) ^ ex PW 

where the values of A and <zk are determined by equations (p"2|). Above the critical temperature in the disordered 
phase the expressions for Tc, Tq and Tq are very similar to the previous ones except for the Lagrange multipliers 
which do not verify Ao = 2DX\ and = 0. In this case the Lagrange multipliers are determined by equations (|^). 
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V. DYNAMICAL SOLUTION 



A. General Method 

The solution to the previous equations is quite involved in the off-equilibrium regime and cannot be exactly cal- 
culated even if some results can be obtained in the asymptotic long time limit. While it is possible to simplify the 
analytical expressions using Laplace transformations the analytical analysis of the dynamical equations appears to 
be tedious. Here we follow a different and more straightforward strategy and numerically investigate the solution to 
the dynamical equations. To understand the nature of the off-equilibrium behavior of the model we have numerically 
integrated the dynamical equations by truncating the hierarchy up to a given finite number of elements. We have 
investigated the cases D =1 and 3. In the first case, there is no phase transition while there is a transition in the 
second one. For the numerical integration of the equations we have considered between 100 and 500 elements of the 
hierarchy for D = 1 and between 5000 and 20000 elements in the condensed phase in D = 3 (where relaxation is 
slower). We note that all figures in this section when plotted in a logarithmic scale are always in base log 10 . In all 
the cases, the integration was performed with an Euler method of second order. Some of this results were tested with 
a fourth order Euler method. 



B. One time quantities 

1. Relaxation of the density 

The first behavior we can study is the relaxation of the density to its equilibrium value. In one dimension we find that 
the relaxation is exponential with time for T different from zero as expected in the disordered phase. For T — 0, iwe 
find different behaviors depending on the initial condition as found in the spherical Sherrington-Kirkpatrick model". 
If the initial condition has a macroscopic projection on the equilibrium state then the relaxation is exponential. If 
the projection on the equilibrium state is zero, then the density of the system relaxes to the equilibrium with a power 
law, t -1 . 

In three dimensions, starting from a non-equilibrium initial condition the system relaxes exponentially fast above 
the critical temperature in the disordered phase. In the condensed phase the relaxation is algebraic with aT 1 power 
law (figure 1). The behavior p£the density in the LP model is very similar to the relaxation of the energy in the 
disordered spherical SK modeO. 
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FIG. 1. Time dependent density in 3D p eq — p (from top to bottom) T = 0.0001, 0.1,1 and 10, with homogeneous initial 
condition. The equilibrium densities at T = 0.0001 is 0.499950, at T = 0.1 is 0.454950, at T = 1. is 0.309017 and at T = 10. is 
0.256246. 








2. Hysteresis effects. 



We have performed some temperature cycling experiments in our system. Starting from a random high temperature 
configuration (for D — 3 we start above the condensation transition temperature) wc let the system equilibrate at this 
temperature and later on we decrease the temperature at a finite rate. As far as we are mainly interested in the slow 
cooling regime, to integrate the equations we change the temperature in the differential equations on constant steps 
of AT. The cooling rate is given by r = 4^, where t* is the time the system spends in a given constant temperature. 
We observe that the system departs from the equilibrium line at a temperature T* which decreases as the cooling 
rate decreases. The inverse of the cooling rate yields an estimate of the relaxation time at that temperature T* . 
Below T* the system fails to relax to the equilibrium and remains at a density lower than the equilibrium value. Once 
we reach zero temperature, we start to increase it at the same rate. We find that the non-equilibrium line crosses 
the equilibrium one which indicates that the system also fails to relax to the equilibrium but now it remains at a 
density higher than the equilibrium value. Finally, the non equilibrium line merges again the equilibrium line at a 
temperature of the same order of T*. This behavior can be observed in 1-D (figure 2) as well as in 3-D (figure 3). In 
these experiments there is no apparent difference between the 3-D and the 1-D case indicating that this is a general 
non-etwilibrium effect unsensitive to the existpace of a phase transition. Similar effects are observed in the 1-D Ising 
modell3 as well as in the Backgammon modeltS. 
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FIG. 2. Cooling experiment in 1-D. 0.5-density versus temperature. We show the equilibrium value for density (continuous 
line) and the value for the density at different cooling rates. These are (from top to bottom) 0.1,0.01 and 0.001. 
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FIG. 3. Cooling experiment in 3-D. 0.5-density versus temperature. We show the equilibrium value for density (continuous 
line) and the value for the density at different cooling rates. These are (from top to bottom) 0.1,0.01 and 0.001. 



3. Heating Experiments 



An striking feature of this system is that although no potential energy has been introduced there appear to be a 
large number of stable crystalline states at zero temperature which have an infinite lifetime. To see the effect of these 
crystalline states on the dynamics we have done some heating experiments. At T — an in 1-D, we have constructed 
different periodic crystalline initial conditions (putting the density every p sites equal to one). As we have shown in 
section (II B), crystalline configurations with p = 2, 3 are metastable states at zero temperature. Numerically we find 
that the configurations where p = 4 and p = 6 decay to the configuration p = 2, and that p = 5 and p = 7 does not 
decay to p = 2 or p = 3, but to solutions with densities between i (p = 2) and ^ (p — 3). These results are simply 
explained noting that, although the restriction o-i = is a global one, apparently the only way to fulfill this restriction 
is imposing that every occupied site is surrounded by empty sites. 

At finite temperature, the metastable states have a finite lifetime and decay to the equilibrium state. However, at 
low temperatures, the lifetime of metastable states can be very large and their effects on the dynamics can then be 
observed. Suppose we prepare the system in such a way that the initial condition belongs to the basin of attraction 
of one of these stationary states. If we perform a heating experiment in which we raise the temperature of the system 
at a finite rate, then we find that the system is trapped for some time near the basin of attraction. In figure 4, we 
show the case where we start from a non metastable state with p — 5. We see that first it decays to the crystal state 
(the first plateau seen in the figure) an later on it relaxes to the equilibrium line. 
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FIG. 4. Heating experiment in 1-D. 0.5-density versus temperature. Starting from a 5 period 1-D crystal initial condition. 
We show the equilibrium value for density (continuous line) and the value for the density at different cooling rates. These are 
(from top to bottom) 0.1,0.01, 0.001 and 0.0001. 



4- Kinetic growth. 



In the LP model, contrarily to the case of the spherical Sherrington-Kirkpatrick modelEj, it is possible to investigate 
kinetic growth phenomena, i.e. how the condensed domains grow as time goes by. This is interesting because allows 
to ascertain how relevant is dimensionality in non equilibrium phenomena. Let us consider the densities in Fourier 
space eq.(|l7|). Using equation we obtain its time evolution, 



ld\mi 

2 dt 



N6(k)p - A (2| j 5(k)| 2 - NS(k) P ) - AX.D |p(k)| 2 7 (fe) + T 



(62) 



Integrating this equation with the values for the Lagrange multipliers obtained from the hierarchy of the density, 
we get the Fourier transformed densities for the infinite system. During the evolution from a non-equilibrium initial 
condition, we see how a peak around k — tt is formed (in 3-D, around k = (71", tt, it)). As temperature decreases the 
peak becomes sharper. For a better understanding of the results, we can appeal to the equilibrium relations. Using 
(p3|), ( p7| ) and ( p8| ) we can calculate an effective temperature for any |p(k)| 2 . Obviously, in equilibrium, all of them 
will have the same effective temperature, the equilibrium one. 

In figure 5, we show the values of the effective temperature for values of k = k (1, 1, 1), with k ranging from to 
7T in the 3-D case (the figure is symmetric around k = n). For times larger than a typical time t* ~ 1 a plateau 
appears in the effective temperature and its value decays very fast to T, the temperature of the heat bath. For values 
of k ranging from zero up to-a given value k maxi all the modes have approximately the same effective temperature 
which is the equilibrium oneEEI. There is a given value of k (let us call it k* ) where the effective temperature becomes 
infinite and above this value of k max the effective temperature is not denned (This means that there does not exist 
an equivalent equilibrium system described by the dynamical densities |p(k)| 2 ). 
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FIG. 5. Effective temperatures versus k in 3-D and at T = 0.1. It is shown for different values of time: 0.01,0.1,1,10,100,1000. 
The vertical line on the right of the plot is k — tt. 



With the value of k max we can also estimate the correlation length using the relation £ 



We find that the 



correlation length diverges as t ? typical of domain growth in ferromagnets with non conserved dynamics (figure 6^-, 
The dynamics in the LP model can be intuitively understood in the framework of the kinetic growth scenarioeil. 
For times less than t* there is no plateau in the effective temperature because the system is nearly empty and the 
system is filling the lattice in a random and uncorrelated way. As soon as the density becomes large enough, any 
change of the density in a site of the lattice is correlated with that of the nearest neighbors and any change in the 
density requires a reorganization of the local densities inside a given domain. So, a critical and finite time t* is 
needed till the first domains appears. At times larger than t* correlated domains start to grow in time. At length 
scales smaller than the growing correlation length £ the system is in local equilibrium (the effective temperature is 
the equilibrium temperature) while it is completely disordered (the effective temperature is infinite) at larger length 
scales. The domain growth scenario is nicely reproduced in the LP model. 




2.0 

log(t) 



FIG. 6. Inverse of the correlation length, l/£, against time in 3-D and at T = 0.1 (dotted line) and 1 (dashed line). The 
straight lines show t -1 ^ 2 behavior. 
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C. Two Times Quantities. 



1. The Correlation and Response Functions. 



In this section, we are going to study the correlation and response functions in equations ( f47| ) and ( pLq ) with x = 0. 
We consider the following normalized correlation function, 



C o (t,t',x = 0)-p 2 eq _ C (t,1f,x = ) - pi 
C (f, f, x = 0) - pl q p - pl q 



C norm (t,t') — — y~ — 5 ~ (63) 



where po is the initial density at time t' and p eq is the equilibrium density corresponding to the system at temperature 
T. In what follows we redefine the time variables t' — t w and consider the behavior of correlation functions for 
different values of t w . 

Our results confirm the simplest Mode Coupling scenario where relaxation proceeds in two steps, the a and f3 
relaxation processes. A comment now seems to be appropriate. Originally Mode Coupling theory was devised ta 
understand the equilibrium dynamics of liquids in the vicinity of the glass transitions. Gotze and collaboratorsO 
have proposed a general mathematical framework for the understanding of equilibrium relaxational processes which 
take place in the vicinity of a bifurcation instability. The model we propose here is a very simplified version of 
the glass scenario where there is no discontinuity of the ergodicity parameter in the bifurcation point (i.e. at the 
condensation phase transition temperature). In fact, above Di — 2 the LP model exhibits a usual second order phase 
transition with classical critical exponents v = 1/(D — 2), rj = (D < 4). The dynamical critical exponent is z = 2 
typical of mean-field theory. This implies that the relaxation time diverges close to the transition temperature T c 
like r ~ (T — T c )~ 2 ' ( D ~ 2 ), In the critical point equilibrium time correlations decay like t~( D ~ 2 )/ 2 , i.e. like £ -1 / 2 for 
D = 3. Among other results these scaling relations led naturally to the \jt decay found in figure 1 for the density 
in all dimensions at the critical point and below it. Above T c the scaling behavior C{t) ~ t~( D ~ 2 '' z f(t/r) holds 
where r is the divergent characteristic time scale and the superposition principle is then valid. We want to stress 
that, contrarily to other scenarios for glassy relaxations (see the discussion by GotzeEj) in this case there is only 
one relaxation process above T c since there is not discontinuity in the ergodicity parameter at T c . But still it is 
quite interesting to investigate the extension of the Mode Coupling scenario to the off-equilibrium dynamics below 
T c . In this case, it is necessary to consider the initial time dependence in the dynamical equations. This implies 
the emergence, among others, of new off equilibrium phenomena like aging, i.e. the dependence of the evolution of 
the system on the initial time state. Also below T c off-equilibrium correlation functions are expected to display the 
characteristic two steps form since the ergodicity parameter (i.e. the Edwards- Anderson parameter) is already finite. 

We will consider t w larger than t* , i.e. the typical time needed to reach a macroscopic density in the lattice starting 
from a nearly empty lattice. For values of t w < t* , the system is quite far from the asymptotic long time regime and 
obviously the two step form is hardly seen. There are qualitative differences between the 1-D and 3-D cases since in 
3-D there is a condensed phase while the system is always disordered in the 1-D case. The off-equilibrium correlation 
function in 3-D decays in two steps: the (3 process or stationary part with a decay to a plateau with a t -1 / 2 and 
the slow a process where the density-density correlation function decays like t~ 2 . In the process the relaxation is 
stationary and depends only on the time difference t — t w while in the slow a process there is aging and the correlation 
functions depend on both t and t w (figure 7). The scaling behavior t/t w is well reproduced in the off-equilibrium 
regime (figure 8). In 1-D the behavior is similar except for the plateau which is absent (only a small inflexion for the 
correlation function at short times is observed). In 3-D the plateau persists over an arbitrarily long time scale which 
grows with t w whereas for 1-D aging is interrupted and disappears when t w is of the order of the finite relaxation 
time. 
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FIG. 7. Correlation function versus t — t w in 3-D at T = 0.1. Different waiting times are shown. From top to bottom 
t w = 10000,1000,300,100,30,10,3,1. 
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FIG. 8. Correlation function versus (t — t w )/t w in 3-D and T=0.1. Different waiting times are shown, 

t w — 1000,300,100,30,10,3,1. For times larger than 1, all of them merge in one curve except for t w = 1 (the lines on the 
left) which is too short. 

The dynamical scenario in the LP model can be inferred from a study of the response function. In any dimension 
the response function decays very fast to zero showing no aging in the asymptotic long time regime. Indeed, in 3D 
we see that G(t,t w ) decays like t~ 3 / 2 for long times and then the integrated response function decays as t^ 1 / 2 . This 
is an indication of the short-time memory of the-system. In the context of glassy dynamics this indicates that the LP 
model has no anomaly in the response functioned. 
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FIG. 9. Response function, G(t,t w ), versus t — t w in 3-D and for T = 0.1. Different waiting times are shown. From top to 
bottom t w = 10,100,1000,10000. 

D. The Fluctuation Dissipation Theorem. 

Once we know the values of the correlation and response functions, we can study the fluctuation-dissipation ratio, 

at™ 

This ratio is equal to 1 in equilibrium (recall that, by definition, a factor T has been absorbed in the response 
function). Roughly, we find that X(t,t w ) is approximately equal to 1 for t — t w < t w , showing that the system is in 
local equilibrium in the stationary regime. For times t — t w larger than t w the X decays to zero very fast. We find 
the same qualitative behavior in 1-D and 3-D. This is expected in the absence of anomaly since the response function 
decays very fast to zero in this regime. If we look more carefully (see figure 10), we find that there are two different 
regimes. For small values of t w (t w < t* ~ 1 when the lattice starts to be filled of particles), the value of X decreases 
monotonically; but for larger values of t w , X increases with t until it reaches a maximum and after decreases very fast. 
Numerically, we can extrapolate that for t w tending to infinity and t — t w <t w (i.e. in the /3-regime) the 1 — X(t, t w ) 
tends to zero as (t — tw) -1 . 

Using the fact that for long times the response functions decays as t~ 3 / 2 and that the X(t,t w ) ~ 1, we can conclude 
that the correlation function decays to the plateau with a power law i~ 1/<2 . This result has been checked directly by 
fitting the decay of the correlation function to the plateau (see below for an estimate of the value of the plateau) for 
large values of t w . 



1G 



0.0 2.0 
l°8(t-U 



4.0 



FIG. 10. 1-X(t, t w ) versus t-t w . Different waiting times are shown. From top to bottom t w =l,3, 10,30, 100,300, 1000, 10000. 
The straight continuous line shows t _ t behavior. 

We can also study the evolution of X(t, t') as a function of C(t, if). Qualitatively we find very similar results in all 
dimensions: in the /3-regime X is close to 1 while in the a-regime it decays very fast to zero. The 1-D and 3-D cases 
are depicted in figures 11 and 12 respectively. Interestingly enough in the 3-D case we also find (figure 12) that all 
the curves (except t w — 1 which is a too short time) cross at C* — C norm {t,t w ) ~ 0.75 (the value of C* increases for 
lower temperatures) corresponding to a value of X* ~ 0.25 . The value C* corresponds to the value of C norm (t w ,t) in 
the plateau in the /3-regime (see figure 7). It is possible to make all different curves (corresponding to different values 
of t w ) collapse in a universal one. We find that the scaling law, 

X(t, t w ) = X((C* - C norm (t, t w )) t a w 4 ) (65) 

for C norm (t w ,t) < 0.75 nicely fits data for all different values of t w (inset in figure 12). We can interpret naturally 
the value of X in the two regimes as a ratio between the temperature of the system and an effective temperature Tf. 
For C norm (t w ,t) > C* we find X ~ 1 and the effective temperature of the system coincides with the temperature 
of the thermal bath. In this /3-regime the system is in local equilibrium. For C norm (t w ,t) < C* we find X ~ 
which means that the effective temperature is infinite. This is a confirmation of the results obtained in section V.B 
interpreted within the kinetic domain growth scenario. In that case it was found that the effective temperature outside 
equilibrated domains was also infinite. Note that the definition of an effective temperature is not free of inconsistencies 
in the most general case and there is not evidence a priori that both effective temperatures (in the a regime) for the 
one time quantities and the effective temperature obtained for the two time quantities coincide. A result of this bfe 
was found in the Backgammon model for the physical interpretation of the violation fluctuation dissipation ratioci 
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FIG. 11. X(t,t w ) versus C n0 rm(t,t w ) in ID at T = 0.0001. Different waiting times are shown. From top to bottom 
t w = 1000,300, 100, 10. 
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FIG. 12. X(t,t w ) versus 1 — C n0 rm(t,t w ) in 3-D at T — 0.1. Different waiting times are shown. From bottom to top 
t w = 1, 3, 10, 30, 100, 300, 1000, 10000. In the inset we show the collapse of the data using the scaling relation eq.(|65|) for various 
values of t w = 30, 100, 300, 1000. 



1. Correlation Between Replicas. 

Up to now we have found qualitatively similar results in the off-equilibrium behavior in 1-D and 3-D. The question 
arises if it is possible to infer from the dynamics the existence of a condensed, growing phase which could distinguish 
the 1-D from the 3-D behavior at finite T. It has been recently suggestecO that it is possible to characterize the 
dynamical behavior in terms of the overlap between two replicas which start in the same initial configuration and 
follow different noise realizations. We have analyzed the quantity, Q(t w ,t) — Q eq where Q eq = (? eq for different values 
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of t w . In a disordered phase we expect that Q(t w , t) — Q eq should decay to zero for long times because the two replicas 
depart one from each other even if they are at the same initial condition at t w . This is the behavior we find in ID 
(figure 13). In the 3-D case in the condensed phase (see figure 14) the situation is different. Now the two replicas 
remember they were in the same configuration at t w and do not depart indefinitely one from each other. The system 
is then constrained to follow something similar to gutters or channels in phase spaceE-1. Intuitively it is not difficult 
to interpret this result in terms of a domain growth process. We already now that the dynamics of the LP model is 
essentially dominated by the growth of correlated domains. In the 1-D case domains appear and disappear in time 
because the system is in the disordered phase. In this case the typical length of domains grow in time (we are at 
low temperatures) but domains can always appear and disappear. The pattern structure of domains is in some sense 
continuously changed in time. In the 3-D case, once domains start to grow, they are not destroyed again. During 
the condensation process the pattern structure of domains is essentially unchanged and it is only rescaled in time. 
Consequently the two-replicas overlap Q(t w ,t) — Q eq does not decay to zero for long times. 
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FIG. 13. Correlation between replicas, Q(t,t m ) — Q eq against t — tw in ID and T = 0.0001. Different waiting times are 
shown. From top to bottom t w = 1,10,100,1000. 
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FIG. 14. Correlation between replicas, Q(t,t w ) — Q eq against total menus waiting time 3-D and T — 0.1. Different waiting 
times are shown. From bottom to top t w = 1,3,10,30,100,300,1000. 

VI. CONCLUSIONS 

We have analyzed in detail the dynamical properties of the LP model originally introduced to understand the 
thermodynamic properties of hard spheres lattice gases in the spherical approximation. From the viewpoint of the 
dynamics this is an interesting model because it is exactly solvable allowing a detailed investigation of the off- 
equilibrium scenario. The model has no built-in disorder and slow dynamics appears as a consequence of the short 
range dynamical constraints present in the system. The dynamics of this modeLshares a large number of features 
with the off-equilibrium dynamics of the spherical Sherrington-Kirkpatrick modeO where dynamics is driven by the 
macroscopic condensation of the system on the disordered ground state. This gives support to the result that disorder 
is not an essential ingredient for the dynamical glassy scenariocm to be valid. 

We have presented a detailed investigation of the dynamical equations of the model by considering the one time 
and two time quantities. One advantage of the model is that_Lt includes short-range effects which are not present in 
other type of mean-field models like the spherical SK modelij. If the initial configuration has density far from its 
equilibrium value then there is a short-time regime where the system is filled very fast in a spatially uncorrelated way. 
The typical time t* for this filling process is of order 1. It is only after that time that slow relaxation starts when the 
system is spatially correlated and needs to reorganize large regions in order to increase its density. 

For values of time larger than t* , we see that the dynamics oL-the model shows striking similarities to the kinetics 
of domain growth of ferromagnets with non conserved dynamicsO. This is expected since the dynamics proposed in 
cq.(^) corresponds to the coarsening dynamics in the so called model B in phase ordering kinetics. In particular we 
have seen that for length scales smaller than a typical correlation length £ (which grows in time) the system is in local 
equilibrium and spatial fluctuations are determined by the equilibrium temperature. Above the correlation length £ 
the system is completely disordered, hence the typical temperature associated to spatial fluctuations is infinite. The 
correlation length £ is an accurate measure of the typical size of the growing domains. It would be interesting to 
investigate the dynamics of the LP model with the conserved dynamics of kinetic growth. In this case, we expect that 
a similar scenario would apply but with different exponents. 

Further support to this domain growth scenario has been obtained from the study of the two time quantities. In 
particular we find that off-equilibrium relaxation proceeds in two well defined steps: a fast j3 relaxation process where 
the fluctuation dissipation relation is obeyed, followed by a slow a process where time translation invariance is lost and 
the fluctuation dissipation ratio is zero. The physical meaning of the a and f3 process is quite appealing in terms of 
domain growth kinetics. The fast /3-process is associated to local rearrangements of densities inside domains whereas 
the slow a-process is associated to the growth of the typical size of these domains. On the other hand, the response 
function shows no aging and decays very fast to zero. This is the scenario of glassy dynamics without anomaly in 
the response function typical of glassy systems with short term memory. A study of the replica-replica overlap has 
also revealed that the growth mechanism is different in the presence or absence of condensation transition. In 1-D 
where there is no condensation phase transition (the system is always in the disordered phase) the domains appear 
and disappear in time even if their typical size increases. In the 3-D case the domains always tend to grow and the 
pattern structure of domains is always maintained. In this last case we find that lim^—^ lim^oo Q(t w , t) is finitec-X 

We want to note that the relaxation of the density as well as correlation functions do not display the phenomena 
of stretching characteristic of glasses. This is due to the oversimplification inherent to the model where only entropy 
barriers are introduced through global constraints (while in the classical hard-spheres model constraints are always 
local). It would be very interesting to introduce some kind of local constraint which could restore the main features 
observed in real glasses. In this direction, it would be quite interesting to extend this research by considering the 
LP model with the dynamics recently proposed Jay Dean where local densities can never become negative, a feature 
which is not considered in the present dynamicsEj. This requires the introduction of noise correlated with the local 
densities, a kind of local constraint. Also it would be interesting if such a dynamics could be exactly solved. 
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